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ABSTRACT 

Using the spectra of a sample of ~28,000 nearby obscured active galaxies from Data Release 7 of the Sloan Digital 
Sky Survey (SDSS), we probe the connection between active galactic nucleus (AGN) activity and star formation 
over a range of radial scales in the host galaxy. We use the extinction-corrected luminosity of the [O m] 5007 A 
line as a proxy of intrinsic AGN power and supermassive black hole (SMBH) accretion rate. The star formation 
rates (SFRs) are taken from the MPA-JHU value-added catalog and are measured through the 3" SDSS aperture. 
We construct matched samples of galaxies covering a range in redshifts. With increasing redshift, the projected 
aperture size encompasses increasing amounts of the host galaxy. This allows us to trace the radial distribution 
of star formation as a function of AGN luminosity. We find that the star formation becomes more centrally 
concentrated with increasing AGN luminosity and Eddington ratio. This implies that such circumnuclear star 
formation is associated with AGN activity, and that it increasingly dominates over omnipresent disk star formation 
at higher AGN luminosities, placing critical constraints on theoretical models that link host galaxy star formation 
and SMBH fueling. We parameterize this relationship and find that the star formation on radial scales <1.7 kpc, 
when including a constant disk component, has a sub-linear dependence on SMBH accretion rate: SFR oc M 036 , 
suggesting that angular momentum transfer through the disk limits accretion efficiency rather than the supply from 
stellar mass loss. 
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1. INTRODUCTION 

Supermassive black holes (SMBHs) are intrinsically linked 
to their parent galaxies. As the black hole grows via accretion in 
the active galactic nucleus (AGN) phase, the bulge of the host 
galaxy builds up from star formation, as evidenced by the corre- 
lation between SMBH mass and velocity dispersion of the host 
galaxy (Kormendy & Richstone 1995; Ferrarese & Merritt 2000; 
Gebhardt et al. 2000; Tremaine et al. 2002; Graham 2012; 
Graham & Scott 2013). Although the evolution of SMBH accre- 
tion and star formation are surely intertwined, the mechanisms 
that couple these processes are not definitively known. Gas can 
be driven to the potential well of the black hole by galaxy merg- 
ers, which compress gas and thereby trigger host galaxy star 
formation (e.g., Sanders et al. 1988; Kauffmann & Haehnelt 
2000; Di Matteo et al. 2005; Springel et al. 2005; Hopkins et al. 
2005). Alternatively, secular mechanisms such as torques from 
galactic bars, bars-within-bars, and spiral arms (Shlosman et al. 

1989; Garcfa-Burillo et al. 2005; Hopkins & Quataert 2010); 
stochastic accretion of cold gas (Hopkins & Hernquist 2006); or 
mass loss from stellar winds (Norman & Scoville 1988; Ciotti 
& Ostriker 1997, 2007; Kauffmann & Heckman 2009) can con- 
tribute significantly to both processes. Feedback is another key 
ingredient in the star formation and AGN connection, where 
AGN winds can blow out circumnuclear gas, suppressing star 
formation and regulating SMBH and bulge growth (Silk & Rees 
1998; Haehnelt et al. 1998; Granato et al. 2004; Di Matteo et al. 
2005; Wyithe & Loeb 2003; Somerville et al. 2008). Observa- 
tional studies of local optically obscured AGN (Type 2 Seyferts, 
Sy2s) have indicated that star formation processes can both 
feed the black hole, via stellar mass loss, and limit accretion 
efficiency through supernova feedback (Wild et al. 2010). 


Seyfert 2 galaxies are prime candidates for studying the 
observational link between SMBH growth and host galaxy star 
formation. As the accretion disk is enshrouded by circumnuclear 
dust obscuration, optical signatures of star formation in the host 
galaxy as well as AGN fueling can be investigated. Previous 
studies of Sy2s have revealed that a significant fraction live in 
host galaxies that are undergoing (or have recently experienced) 
significant amounts of star formation (Kauffmann et al. 2003a; 
Cid Fernandes et al. 2004; Netzer 2009). Using a sample of over 
20,000 AGNs from the Sloan Digital Sky Survey (SDSS) Data 
Release One (Stoughton et al. 2002), Kauffmann et al. (2003a) 
demonstrated that the trend between star formation activity 
and AGN accretion was not confined to the nuclear regions of 
the host galaxies for luminous AGN. Subsequently, Kauffmann 
et al. (2007) used a combination of SDSS spectra and GALEX 
UV images to show that the link between star formation and 
black hole growth was much more direct for star formation in the 
central few kpc than for star formation in the global galactic disk. 
In agreement with this result, Diamond-Stanic & Rieke (2012) 
studied 74 galaxies from the revised Shapley-Ames Seyfert 
sample (Shapley & Ames 1932; Maiolino & Rieke 1995; Ho 
et al. 1997) and found that this connection is strongest on sub- 
kpc scales, with a much weaker correlation when considering 
larger scale (i.e., > 1 kpc) star formation. These place important 
constraints on the mechanisms by which star formation and 
black hole growth are linked. 

However, these studies suffer from non-uniform methodology 
for measuring star formation rates (SFRs) and /or small sample 
sizes. To address these limitations, in this work we test the 
relationship between star formation and SMBH accretion on 
various galactic scales using a sample of ~28,000 type 2 AGN 
and composite galaxies culled from SDSS. All of the quantities 
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are measured in a uniform way through the 3" spectroscopic 
fiber. These quantities parameterize the AGN luminosity, SMBH 
accretion rate, and SFR. We test whether the connection between 
AGN fueling and starburst activity varies as the projected size of 
the SDSS aperture sequentially covers greater parts of the host 
galaxy, and whether this relationship is a function of intrinsic 
AGN luminosity. 

2. DATA ANALYSIS 

Our data were drawn from the Main Galaxy Sample of 
SDSS Data Release 7 (Strauss et al. 2002; Abazajian et al. 
2009) where we utilized the MPA-JHU value added catalogs. 4 
We isolated sources categorized as “AGN” and “composite” 
according to the diagnostic Baldwin-Phillips-Terlevich (BPT) 
diagram (Baldwin et al. 1981; Kewley et al. 2001; Kauffmann 
et al. 2003a); we note that these classifications were made 
automatically by the MPA-JHU pipeline. We exclude galaxies 
classified as LINERs. We imposed several quality control cuts: 
the Ha, H/3, [Oiii] 5007 A, and [Nil] 6584 A emission lines 
had to be significant above the 5er level and the Ha equivalent 
width had to exceed 3 A to remove galaxies where the ionizing 
radiation may be due to low-mass evolved stars and not nuclear 
accretion (Cid Fernandes et al. 2011). We also restrict the 
sample to the Petrosian magnitude range 15.5 < m r < 17.8 (the 
limiting magnitude of the Main Galaxy Sample) and remove 
outliers in BPT parameter space, focusing between —1.0 < 
log([N n]/Ha) < 1.0 and -1.0 < log([0 iii]/H/3) < 2. The 
objects in our sample span 0.01 < z < 0.32. As shown by 
Kauffmann et al. (2003a) and Reichard et al. (2009), these 
AGN host galaxies are typically structurally normal galaxies of 
intermediate Hubble type (average SDSS concentration index 
of 2.61 ±0.33). 

2.1. AGN Activity 

We use the luminosity of the [Oiii] line (L[Om|) as a proxy 
for intrinsic AGN luminosity since it is primarily ionized by 
accretion disk photons. Since [Oiii] forms in the AGN nar- 
row line region, it is not affected strongly by dust obscuration 
from the circumnuclear torus, providing a good estimate of the 
intrinsic AGN luminosity in SDSS galaxies (Kauffmann et al. 
2003a; Heckman et al. 2004; LaMassa et al. 2010). For compos- 
ite galaxies, i.e., those systems between the Kauffmann et al. 
(2003a) empirical boundary between star-forming and active 
galaxies and the Kewley et al. (2001) theoretical maximum 
starburst line on the BPT diagram, we estimated the AGN con- 
tribution to the [Oiii] emission following the prescription of 
Wild et al. (2010). This method estimates the star formation 
contributions to L[Om] based on the location of the source in 
the BPT diagram. We only keep composite galaxies if the AGN 
component of the [O in] line is at least 30% of the total [O hi] 
emission, removing very weak AGN. In total, we have 27,952 
galaxies, of which 15,508 are composites and 12,444 are Sy2s; 
they are shown on the BPT diagram in Figure 1. Finally, we 
correct Z,[Om] for dust obscuration using the observed Balmer 
decrement (Ha/H/3), assuming an intrinsic Ha/H/i value of 

3.1, and the standard R(V) = A(V)/E(B — V) — 3.1 extinction 
curve (Cardelli et al. 1989), obtaining L[om],corr- 

The SMBH mass (Mbh) goes as the host galaxy stellar 
velocity dispersion (er) to the fourth power (Tremaine et al. 
2002). To parameterize the SMBH accretion rate, we therefore 


4 http://www.sdss3.org/dr9/algorithms/galaxy_mpa jhu.php 
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Figure 1. BPT diagram for galaxies used in our analysis. The lower dashed line 
represents the empirical boundary between star-forming galaxies and active 
systems from Kauffmann et al. (2003a) while the upper dotted line delineates 
the maximum theoretical starburst line from Kewley et al. (2001). There are 
15,508 composite galaxies (red data points, between the two curves) and 12,444 
Sy2 galaxies (blue data points). These composite systems represent the subset 
of composites in which the AGN contribution to [O hi] emission is estimated at 
^30% of the total. 

(A color version of this figure is available in the online journal.) 

define a pseudo-Eddington parameter, L[o nl ] C0rr /cr 4 , for the 
galaxies where a velocity dispersion was measured. In the 
following, we search for trends between AGN luminosity 
(£[Oiii],corr) and SFR, and between the Eddington parameter 
(E[Oin],corr/o’ 4 ) and specific SFR (sSFR, the SFR normalized by 
galaxy mass). 

2.2. Star Formation 

We focus on the SFR and sSFR measured through the 
3" SDSS spectroscopic fiber (SFRfiber, sSFR^ber), calculated 
as outlined in Brinchmann et al. (2004). As nuclear activity 
contributes to Ha luminosity, which is used to calibrate the SFR 
in quiescent galaxies, these SFRs are instead calculated using the 
break in the optical spectrum at 4000 A (D„(4000); Balogh et al. 
1999), which is a measure of the age of the host galaxy stellar 
population: a shallower break is associated with a younger stellar 
population. Using the relationship between Z)„(4000) and the 
Ha calibrated sSFR^ber for star-forming galaxies, Brinchmann 
et al. (2004) constructed probability distribution functions from 
which sSFRfiber were estimated for active galaxies; SFRgber is 
then calculated by multiplying sSFR^ber by the galaxy mass 
within the fiber. In this analysis, we use the value-added 
catalog provided median estimates of the SFR| lhcr and sSFR| lhcr 
probability distribution functions. 

2.3. z/Zmax as a Distance Proxy 

Our goal is to determine the relationship between black hole 
growth and star formation over a range of physical scales inside 
the galaxy. To ensure that we are comparing host galaxies with 
similar properties at all distances, isolating the physical size 
of the projected SDSS fiber aperture as the varying parameter, 
we introduce the z/z m ax parameter (Kauffmann et al. 2003a, 
2003b), where z max is the redshift at which the galaxy falls out 
of the apparent magnitude limit of the sample ( m r — 17.8). 
Absent significant cosmic evolution (see below), the properties 
of a population of galaxies selected from SDSS will be invariant 
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Figure 2. (a) SFRfiber a s a function of AGN luminosity and (d) sSFRfiber as a function of Eddington parameter in increasing z/z m ax bins, (b, e) associated dispersions, 
and (c, f) legends. The legends give the median redshift of the galaxies in each bin: with increasing redshift, the projected physical aperture size grows, encompassing 
greater amounts of the host galaxy. The relationship between AGN luminosity and star formation shows a dependence on the amount of the host galaxy that is sampled 
at lower AGN luminosities, but not at higher ones. The slope of the relation between sSFR and Eddington parameter becomes flatter as the aperture size grows, 
meaning that the star formation in more luminous AGN is centrally concentrated. 

(A color version of this figure is available in the online journal.) 


with z/zmax- Normalizing the redshift of each galaxy by z max 
corrects for the effect that more luminous and massive galaxies 
are visible out to greater distances. To calculate z max , we first 
calculate the absolute r magnitude (M, ) using 

M r = m r — DM(z) — k(z), (1) 

where DM(z) is the distance modulus and k(z) is the 
^-correction. We adopted a cosmology of Hq — 70 km s 1 Mpc - 1 , 
Qm = 0.27, and A =0.73. 

We determined k-corrections for each galaxy using kcorrect 
version 4.2 (Blanton & Roweis 2007). Using M R , we then found 
the redshift, Zmax, at which each galaxy would have an apparent 
magnitude of 17.8. As indicated in the legends of Figure 2, 
z/Zmax is proportional to redshift. 

2.4. The Effects of Cosmic Evolution and Sample Completeness 

With our approach of using galaxies over a range of redshifts 
to study the radial distribution of star formation in AGN, there 
are two issues that must be addressed. First, since we are 
confining our study to the relatively local universe, we do not 
expect the number density of AGN to evolve with redshift. Our 
sample binned in z /z m ax only spans a range in median redshift 
from 0.06 to 0. 13. Over this range, the increase in the rate of star 
formation and black hole growth per comoving volume element 
(e.g., Marconi et al. 2004) is only a factor of ~1.2 (0.08 dex). 
Second, low-luminosity AGN will be increasingly difficult to 
recognize in more distant galaxies since the diagnostic emission 
lines will become fainter and more diluted by emission due to 
host galaxy star formation. As we noted in LaMassa et al. (2012), 
mis-identification of AGN as star-forming galaxies in SDSS at 


Table 1 

Median V / V ma x for AGN Luminosity and Eddington Parameter Bins 


(L[Om],corr) 

(erg s -1 ) 

WVnax 

(^[ Om],corr/ or ^) 
(erg s 3 km -4 ) 

V/Vmax 

38.5-39.0 

0.25 

30.0-30.5 

0.36 

39.0-39.5 

0.33 

30.5-31.0 

0.32 

39.5-40.0 

0.38 

31.0-31.5 

0.37 

40.0-40.5 

0.44 

31.5-32.0 

0.43 

40.5-41.0 

0.49 

32.0-32.5 

0.46 

41.0-41.5 

0.51 

32.5-33.0 

0.50 

41.5-42.0 

0.55 

33.0-33.5 

0.53 

^42.0 

0.60 

^33.5 

0.59 


Z < 0.3 is rare, and so we do not expect our results to be biased 
by missing a large population of host-galaxy obscured AGN. 

To address these two issues, we performed a V/V max test 
for our sample, where V is the cosmic volume probed by an 
individual galaxy sample and Vmax is the maximum volume out 
to which that galaxy would be observable before dropping below 
the sensitivity limit of the survey. Values for the median of a 
V / V max distribution of significantly less than (greater than) 0.5 
signal incompleteness (cosmic evolution). In Table 1, we list 
the median V / V max for various AGN luminosity and Eddington 
parameter bins. In the following, we work within the range 
where V/V max = 0.5 ± 0.1, i.e., L [0 m],corr > 10 40 erg s _1 
and L[Om\,con / ) 3 x 10 31 erg s 3 km -4 . We also find that 
the median estimate of the total stellar mass and the velocity 
dispersion does not evolve with z/z m ax in these AGN luminosity 
and Eddington parameter bins. Thus, we are comparing similar 
galaxies in aggregate over increasing z/z max values. 


3 


The Astrophysical Journal Letters, 765:L33 (6pp), 2013 March 10 


LaMassa et al. 



1.5 : 
1.0 ; 
0.5 ; 
o 0.0 ; 
-0.5 ; 
- 1.0 : 


-1.5 

40.0 


(a) 


40.5 41.0 41.5 

log ( L [ow],corr erg/s) 

(b) 


42.0 


42.5 



(c) 


(d) 


Figure 3. (a) Log (SFRgberl as a function of z/z m ax in increasing L[om],coir bins and (b) the best-fit slope in each bin. (c) Log (sSFRfiber) as a function of z/Zmax in 
increasing L[om],corr/o ' 4 bins and (d) the best-fit slope in each bin. At higher AGN luminosities and accretion rates, galactic star formation becomes more centrally 
concentrated. The lines in (a) correspond to the best-fit parameters from SFRfjber = SFR nuc i ea i- + S LRdi sk - where SFR nuc i ear = a (Mq yr 1 ) x(L[om],corr/10 42 erg s ~* f’ 
and SFRdisk is SI Rohcr associated with the lowest AGN luminosity bin (i.e., 40.0 < log(L[om],corr < 40.5 dex). From a reduced y 2 minimization, we find a = 0.44 ± 
0.02 and p = 0.36 ± 0.04. 

(A color version of this figure is available in the online journal.) 


Adopting the bolometric correction for the extinction- 
corrected [O in] luminosity from Kauffmann & Heckman (2009) 
and the conversion from a to M BH from Tremaine et al. (2002), 
our sample covers ranges from roughly 10 92 to 10 12,2 L Q in 
bolometric luminosity and from roughly 10~ 2 8 to 1 in L/LEdd- 
Over these broad ranges, we are reasonably complete and the 
effects of cosmic evolution are modest. 

3. RESULTS 

Separating galaxies into bins of z/z m ax, we plot SFRfiber as 
a function of AGN luminosity (L[Om],corr) and sSFRfiber as a 
function of the Eddington parameter (Lipmi.corr/T 3 4 ), and the 
associated dispersions, in Figure 2. Within each z/z m ax bin, 
we further bin in AGN luminosity and Eddington parameter, 
including at least 20 sources in each bin. We then calculate the 
median SFR and sSFR in these bins. The dispersion ( disp ) is 
the associated la error for a Gaussian distribution, i.e., half 
of the difference between the 83rd and 17th percentiles of the 
SFRfiber and sSFRfiber distributions in each bin. 

As z / Z m ax increases, the 3" spectroscopic fiber covers larger 
physical scales of the host galaxy. The median galaxy red- 
shift ranges from ~0.06 to ~0.13 between the lowest and 
highest .c/z max bin, corresponding to physical aperture radii 


of 1.7-3. 5 kpc. As can be seen in Figure 2, the SFR (sSFR) 
increases slowly but systematically with the AGN luminosity 
(Eddington ratio), within a given z/z max bin. It is also clear that 
there are systematic differences between these relations as a 
function of SDSS aperture size (z/zmax)- At low to moderate 
AGN luminosities, there is a trend for the SFR to grow for a 
given AGN luminosity as the projected aperture size increases, 
although this trend is not clear at higher AGN luminosities. 
Moreover, the slope of the relationship between sSFR and the 
Eddington-ratio proxy flattens as z/z m ax increases. 

To see these trends with distance more clearly, in the left pan- 
els of Figure 3, we plot the SFR (top panel) and sSFR (bottom 
panel) as functions of z /z max , and bin the data by the AGN lumi- 
nosity (top) and Eddington ratio (bottom). These plots demon- 
strate that the radial variations of SFR and sSFR are dependent 
on the AGN luminosity and Eddington ratio, respectively. To 
visualize and quantify how the radial distribution of SFR and 
sSFR varies as a function of increasing AGN luminosity and Ed- 
dington ratio, we have calculated the slope, a, of the correlation 
between 

(S)SFR (X (z / Zmax) 0 (2) 

in each AGN luminosity and Eddington parameter bin. For the 
linear regression fit, we weighted the data points by the error in 
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each bin, taken as disp/V/V, where disp is the dispersion in the 
distribution within the bin (half the difference between the 17th 
and 83rd percentiles) and N is the number of sources per bin. 

The results in the right panels of Figure 3 show the same 
systematic trend. As the AGN luminosity increases, the increase 
in star formation inside the SDSS aperture becomes smaller with 
greater projected aperture size. At lower AGN luminosities, the 
relationship between SFR and increasing projected size of the 
host galaxy is roughly linear, showing that star formation is 
distributed over the whole range in radial scales that we probe 
in these host galaxies. However, at higher AGN luminosities, 
star formation is predominantly circumnuclear: we do not see 
evidence that the SFR increases significantly as the size of the 
region we observe increases from 1.7 to 3.5 kpc (the physical 
size of the fiber radius at the median redshifts from lowest 
to highest z/z m ax bins). Similarly, at high Eddington rates, the 
negative slope between sSFR and projected aperture size further 
suggests dominant circumnuclear star formation: as the size of 
the observed region grows, stellar mass is added, but is not 
compensated by an increase in star formation. 

4. DISCUSSION 

These results indicate that circumnuclear star formation 
associated with AGN activity (SFR nuc ] eus ) only dominates over 
omnipresent disk star formation (SFRdisk) when the AGN is very 
luminous. During such a stage, host galaxy star formation will 
appear to be centrally concentrated. Let us assume 

SFRfiber = SFRnucleus + SFR d i s k- (3) 

The star formation in the nucleus is associated with AGN 
fueling, parameterized by 

SFRnucleus = a{M e yr _1 ) x (L [Onl ] >corr /10 42 erg s -1 /. (4) 

We use SFR fiber from the lowest AGN luminosity bin, 40.0 < 
log (L[Oiii],coit) < 40.5 dex, as the omnipresent SFRdisk- We 
solve Equation (3) for all of the log (Zqomi.corr) > 40.5 dex bins 
in Figure 3(a) simultaneously by minimizing / 2 . Here, we take 
the error on each bin to be disp/V7v, where disp is half the 
difference between the 83rd and 17th percentiles and A refers 
to the number of sources in the bin. We find a — 0.44 ± 0.02 
and p = 0.36 ± 0.04, and plot these fits in Figure 3(a); the 
errors represent the 95% confidence interval, where Ax 2 = 3.84. 
These results indicate that star formation associated with SMBH 
fueling rises sub-linearly with AGN luminosity. The fit could 
be improved with a more complex model in which both the disk 
and circumnuclear SFRs are allowed to have separate functional 
dependencies on the AGN luminosity, but such complexity is 
not really warranted. Our simple toy model is only meant to 
show that the data are consistent with a stronger link between 
the AGN luminosity and the amount of central star formation 
than with larger scale star formation. 

What can be the physical link between these two processes? 
Circumnuclear star-formation in luminous AGN can originate 
from major mergers, as simulations have demonstrated that 
torques produced by galaxy-galaxy interactions fuel SMBHs 
while triggering star formation in the galactic center (Di Matteo 
et al. 2005; Springel et al. 2005; Hopkins et al. 2005). In- 
deed, Treister et al. (2012) have demonstrated that merg- 
ers are responsible for triggering the most luminous AGN 
(Lboi > 10 45 erg s -1 ). However, the majority of the highest 
luminosity AGN in our sample falls below this bolometric lu- 
minosity threshold. Additionally, the lopsidedness of SDSS Sy2 
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host galaxies within the Eddington range we probe (~ 10~ 3 to 1) 
does not drastically increase, suggesting that mergers are not the 
primary mechanism for coupling AGN and star formation ac- 
tivity in this sample (Reichard et al. 2009). Alternatively, Ciotti 
& Ostriker (2007) have shown that recycled gas from dying 
stars near the galactic centers of elliptical galaxies becomes ra- 
diatively unstable, causing circumnuclear star formation while 
also feeding the SMBH. As instabilities caused by bars, bars- 
in-bars, spiral arms, etc. produce galactic-scale star formation 
while funneling gas to the SMBH (Hopkins & Quataert 2010), 
stellar mass loss seems a more likely mechanism for coupling 
AGN activity and circumnuclear star formation in moderate lu- 
minosity Sy2 galaxies. 

5. CONCLUSIONS 

We find that the AGN-star formation connection is largely 
driven by circumnuclear starburst activity in local Type 2 AGN, 
consistent with Kauffmann et al. (2007) and Diamond-Stanic & 
Rieke (2012). 

The size- scale of star formation in local active galaxies 
appears to depend on AGN luminosity: it becomes increasingly 
compact in more luminous AGN. More specifically, the lack of 
an increase in the amount of SFR as the size of the region probed 
by the SDSS spectra increases implies that star formation is 
concentrated in radii below <1.7 kpc in these luminous systems. 
This constrains the mechanism that couples AGN activity and 
host galaxy star formation, favoring theories that link these 
two processes at circumnuclear rather than galactic scales (i.e., 
stellar mass loss rather than instabilities caused from bars, etc.). 

Our results are consistent with the hypothesis that circum- 
nuclear star formation is associated with AGN activity, and 
thus dominates over omnipresent galactic disk star formation 
when the AGN becomes luminous. The dependence of the SFR 
on SMBH fueling is sub-linear: on radial scales <1.7 kpc, 
SFR oc M 0 - 36 . This sub-linearity may indicate that accretion 
is a “demand” rather than “supply” driven process. Angular 
momentum transfer through the disk may be relatively ineffi- 
cient, limiting the accretion rate, regardless of the amount of fuel 
provided via stellar mass loss. The remaining material could po- 
tentially be ejected in winds and/or jets, or perhaps accumulate 
in the obscuring “torus.” 

We thank the referee for a careful reading of this manuscript. 
Funding for the SDSS and SDSS-II has been provided by 
the Alfred R Sloan Foundation, the Participating Institutions, 
the National Science Foundation, the U.S. Department of 
Energy, the National Aeronautics and Space Administration, 
the Japanese Monbukagakusho, the Max Planck Society, and 
the Higher Education Funding Council for England. 

REFERENCES 

Abazajian, K. N., Adelman-McCarthy, J. K., Agiieros, M. A., et al. 2009, ApJS, 

182, 543 

Baldwin, J. A., Phillips, M. M„ & Terlevich, R. 1981, PASP, 93, 5 
Balogh, M. L., Morris, S. L., Yee, H. K. C., Carlberg, R. G., & Ellingson, E. 

1999, ApJ, 527, 54 

Blanton, M. R„ & Roweis, S. 2007, AJ, 133, 734 

Brinchmann, J.. Chariot, S., White, S. D. M., et al. 2004, MNRAS, 351, 

1151 

Cardelli, J. A., Clayton, G. C„ & Mathis, J. S. 1989, ApJ, 345, 245 

Cid Fernandes, R., Gu, Q., Melnick, J., et al. 2004, MNRAS, 355, 273 

Cid Fernandes, R„ Stasinska, G., Mateus, A., & Vale Asari, N. 2011, MNRAS, 

413, 1687 

Ciotti, L„ & Ostriker, J. P. 1997, ApJL, 487, L105 


5 


The Astrophysical Journal Letters, 765:L33 (6pp), 2013 March 10 


LaMassa et al. 


Ciotti, L., & Ostriker, J. P. 2007, ApJ, 665, 1038 
Diamond-Stanic, A. M., & Rieke, G. H. 2012, ApJ, 746, 16 
Di Matteo, T., Springel, V., & Hemquist, L. 2005, Natur, 433, 604 
Ferrarese, L., & Merritt, D. 2000, ApJL, 539, L9 

Garcla-Burillo, S., Combes, F., Schinnerer, E., Boone, F., & Hunt, L. K. 
2005, A&A, 441, 1011 

Gebhardt, K., Bender, R., Bower, G., et al. 2000, ApJL, 539, L13 

Graham, A. W. 2012, ApJ, 746, 113 

Graham, A. W„ & Scott, N. 2013, ApJ, 764, 151 

Granato, G. L., De Zotti, G., Silva, L., Bressan, A., & Danese, L. 2004, ApJ, 
600, 580 

Haehnelt, M. G., Natarajan, P., & Rees, M. J. 1998, MNRAS, 300, 817 
Haring, N„ & Rix, H.-W. 2004, ApJL, 604, L89 

Heckman, T. M., Kauffmann, G., Brinchmann, J., et al. 2004, ApJ, 613, 109 
Ho, L. C„ Filippenko, A. V., & Sargent, W. L. W. 1997, ApJS, 1 12, 315 
Hopkins, P. F., & Hemquist, L. 2006, ApJS, 166, 1 
Hopkins, P. F., Hemquist, L., Cox, T. J., et al. 2005, ApJ, 630, 705 
Hopkins, P. F„ & Quataert, E. 2010, MNRAS, 407, 1529 
Kauffmann, G., & Haehnelt, M. 2000, MNRAS, 31 1, 576 
Kauffmann, G., & Heckman, T. 2009, MNRAS, 397, 135 
Kauffmann, G., Heckman, T. M., Budavri, T., et al. 2007, ApJS, 173, 357 
Kauffmann, G., Heckman, T. M., Tremonti, C., et al. 2003a, MNRAS, 
346, 1055 

Kauffmann, G., Heckman, T. M„ White, S. D. M., et al. 2003b, MNRAS, 
341, 33 


Kewley, L. J.. Dopita, M. A., Sutherland, R. S., Heisler, C. A., & Trevena, J. 
2001, ApJ, 556, 121 

Kormendy, J.. & Richstone, D. 1995, ARA&A, 33, 581 

LaMassa, S. M., Heckman, T. M., Ptak, A., et al. 2010, ApJ, 720, 786 

LaMassa, S. M., Heckman, T. M., Ptak, A., et al. 2012, ApJ, 758, 1 

Magorrian, J., Tremaine, S., Richstone, D., et al. 1998, AJ, 1 15, 2285 

Maiolino, R., & Rieke, G. H. 1995, ApJ, 454, 95 

Marconi, A., Risaliti, G„ Gilli, R„ et al. 2004, MNRAS, 351, 169 

Netzer, H. 2009, MNRAS, 399, 1907 

Norman, C„ & Scoville, N. 1988, ApJ, 332, 124 

Reichard, T. A., Heckman, T. M., Rudnick, G., et al. 2009, ApJ, 691, 1005 

Sanders, D. B„ Soifer, B. T„ Elias, J. H„ et al. 1988, ApJ, 325, 74 

Shapley, H., & Ames, A. 1932, AnHar, 88, 41 

Shlosman, I., Frank, J., & Begelman, M. C. 1989, Natur, 338, 45 

Silk, J., & Rees, M. J. 1998, A&A, 33 1 , L 1 

Somerville, R. S., Hopkins, P. F., Cox, T. J., Robertson, B. E., & Hemquist, L. 
2008, MNRAS, 391,481 

Springel, V., Di Matteo, T., & Hemquist, L. 2005, MNRAS, 361, 776 
Stoughton, C., Lupton, R. H., Bemardi, M., et al. 2002, AJ, 123, 485 
Strauss, M. A., Weinberg, D. H„ Lupton, R. H„ et al. 2002, AJ, 124, 1810 
Treister, E., Schawinski, K., Urry, C. M., & Simmons, B. D. 2012, ApJL, 
758, L39 

Tremaine, S., Gebhardt, K.. Bender, R., et al. 2002, ApJ, 574, 740 
Wild, V., Heckman, T„ & Chariot, S. 2010, MNRAS, 405, 933 
Wyithe, J. S. B., & Loeb, A. 2003, ApJ, 595, 614 


6 


